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Abstract 

Numerical studies of gravitational collapse in isotropic coordinates have recently shown 
an interesting connection between the gravitational Lagrangian and black hole thermo- 
dynamics. A study of the actual spacetime was not the main focus of this work and in 
particular, the rich and interesting structure of the interior has not been invetigated in 
much detail and remains largely unknown. We elucidate its features by performing a nu- 
merical study of the spacetime in isotropic coordinates during gravitational collapse of a 
massless scalar field. The most salient feature to emerge is the formation of a thin shell of 
matter just inside the apparent horizon. The energy density and Ricci scalar peak at the 
shell and there is a jump discontinuity in the extrinsic curvature across the apparent hori- 
zon, the hallmark that a thin shell is present in its vicinity. At late stages of the collapse, 
the spacetime consists of two vacuum regions separated by the thin shell. The interior is 
described by an interesting collapsing isotropic universe. It tends towards a vacuum (never 
reaches a perfect vacuum) and there is a slight inhomogeneity in the interior that plays a 
crucial role in the collapse process as the areal radius tends to zero. The spacetime evolves 
towards a curvature (physical) singularity in the interior, both a Weyl and Ricci singularity. 
In the exterior, our numerical results match closely the analytical form of the Schwarzschild 
metric in isotropic coordinates, providing a strong test of our numerical code. 
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1 Introduction 

Isotropic coordinates have recently been used to study numerically black hole thermodynamics 
during gravitational collapse [H [2] and to study numerically the evolution of metric perturba- 
tions during reheating [3]. These coordinates are convenient for numerical studies of gravita- 
tional collapse because the metric functions are not singular anywhere [HE]. In particular, the 
functions are continuous and finite across the horizon and the entire spacetime is covered with 
one coordinate patch. The metric is in a 3 + 1 decomposition form allowing for a Hamiltonian 
or Lagrangian formulation [3]. Moreover, the coordinate t represents the time measured by an 
asymptotic observer at rest and this feature allows one to relate the gravitational Lagrangian 
to a thermodynamic potential making isotropic coordinates ideal for black hole thermodynamic 
studies pQ. In particular, it was shown numerically that in isotropic coordinates the negative 
of the gravitational Lagrangian approaches the free energy of a black hole at late stages of 
collapse [HE]. 

In a 3 + 1 decomposition [4] a metric which is spherically symmetric takes on the following 
general form: 



There are altogether four metric functions of r and t: the lapse N(r,t), the shift /3(r, t) and 
two functions A(r, t) and B(r,t) for the spatial three-metric. We can use gauge freedom, the 
twofold ambiguity associated with choosing two different coordinates, to reduce the number 
to two metric functions. The isotropic coordinate system is based on the choice /3 = and 
"isotropic gauge" A 2 = B 2 = ip 4 for the spatial metric so that (pQ) reduces to 



ip is called the conformal factor. The spherically symmetric vacuum solution in isotropic 
coordinates is given analytically by [5] 



The event horizon is located at r = GM/2 and the region r>GM/2 covers the exterior static 
region of the Schwarzschild black hole. However, the region r < GM/2 does not cover the 
interior of the Schwarzschild black hole; it covers the static exterior region a second time so 
that metric (J3J) is a double covering of the Schwarzschild exterior. In standard coordinates, 
the Schwarzschild metric has the form ds 2 = -A(R) dt 2 + B(R) dR 2 + R 2 dQ, 2 where A(R) = 
1 - 2GM/R and B(R) = A{R)~ 1 . The coordinate transformation R = r (1 + ^f) 2 yields the 
metric in the isotropic form ([3]). We see that r = GM/2 corresponds to the event horizon 
R = 2GM. The regions r > GM/2 and r < GM/2 both correspond to R > 2GM, the exterior 



ds 2 = -(N 2 - A 2 (3 2 ) dt 2 + 2 A 2 fi drdt + A 2 dr 2 + B 2 r 2 dtt 2 . 



(1) 




(2) 
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Schwarzschild region. The interior region R < 2GM of the Schwarzschild black hole is not 
covered by metric ([3j) . This is clear when we consider the Killing vectors. In the interior 
region of the Schwarzschild black hole, all Killing vectors are spacelike [6j and this implies 
that the interior is nonstationary. In contrast, the metric Q has a timelike hypersurface- 
orthogonal Killing vector in the region r < GM/2 and is therefore static in that region. To 
explore the interior in isotropic coordinates, the metric must be of the time-dependent form 
([2]). During collapse, the metric ([2]) is not singular at the horizon: the functions N and ip 
are continuous across the apparent horizon and finite everywhere. Moreover, the coordinate t 
remains timelike in the interior. Note that the coordinate t in isotropic coordinates represents 
the time measured by an asymptotic observer at rest. Such a time coordinate has not been 
used in most previous numerical studies of the interior during spherically symmetric collapse. 
This includes scalar field collapse in both Painleve-Gullstrand (PG) coordinates [7] and in 
coordinates used for studying the evolution of trapped-surfaces [8] , Einstein- Yang Mills collapse 
using maximal slicing and radial spatial coordinates [9] and charged scalar field collapse in 
"Eddington-Finkelstein" |10j or double- null coordinates 

Numerical studies of gravitational collapse in isotropic coordinates have largely focused on the 
thermodynamics [HE] and to date, a separate analysis of the interior spacetime is lacking. The 
interior in isotropic coordinates has a rich structure and we elucidate it by studying numerically 
the entire spacetime during the gravitational collapse of a massless scalar field to a black hole. 
The most salient feature is the emergence of a thin shell of matter at a radius just below r^, the 
radial location of the apparent horizon (where the lapse function N(r,t) is zero at late times). 
The energy density peaks just inside the apparent horizon (and hence inside the event horizon) 
and the extrinsic curvature undergoes a jump discontinuity across the horizon, the hallmark 
that a thin shell has formed in its vicinity [3]. By tracking curvature scalars, the spacetime can 
be seen to be evolving towards a curvature singularity in the interior (both a Weyl and Ricci 
singularity). The region outside the horizon r > evolves towards the metric (j3J), the static 
vacuum Schwarzschild exterior in isotropic coordinates. The interior tends towards a vacuum 
and corresponds to a collapsing isotropic universe where the trace of the extrinsic curvature is 
spatially constant and increases rapidly with time. The interior is not exactly a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe because it is slightly inhomogeneous. The small 
inhomogeneity plays a significant role in the collapse process as ip — > 0. The areal radius is 
given by R = ip 2 r where r is the radial coordinate in isotropic coordinates. In the interior, tp 
ranges from to a value of 2 at r = GM/2 so that the areal radius ranges from to 2M in the 
interior; the entire spacetime, both interior and exterior, is covered in isotropic coordinates. 
In the work of [U |2], both the exterior and interior region make a contribution to the free 
energy and the entire spacetime is needed to obtain numerically the thermodynamics for the 
Schwarzschild black hole. 

Radial null geodesies in isotropic coordinates have dr/dt equal to zero at the N = two-surface 
(where r and t are the radial and time coordinate in isotropic coordinates). Matter accumulates 
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in a thin shell near the N = surface, which at late times is situated at a finite (non-zero) 
radius r^, the radius of the apparent horizon. The thin shell forms just inside the apparent 
horizon and is seen by a static observer in the interior at r = (note that a static observer 
at r = remains static at all times). The proper radial distance between r = and r = 
tends to zero with time so that the thin shell is "falling" towards r = even though it remains 
at a radial location near r^. The proper radial distance is given by J^ h ip 2 {r,t) dr which is 
approximately equal to ip 2 (t)r) l (since tp is approximately homogeneous in the interior). This 
tends towards zero as ip approaches 0. 

Analytical studies have shown that when matter is thrown into a Schwarzschild black hole, a 
static asymptotic observer sees the matter form into thin "pancakes" at the horizon [12j . The 
formation of thin layers of matter at the horizon is not seen by freely falling observers (called 
by the acronym Frefos in [E]). In isotropic coordinates, the time t is the time measured by a 
static asymptotic observer. The formation of a thin shell at a radius near the horizon is the 
dynamical realization of the thin "pancake" scenario discussed in |12| . As far as we know, our 
work is the first time that a numerical collapse study has explicitly shown and discussed this. 
Moreover, our numerical and hence dynamical results yield new insights on this phenomena. 
An analytical test-particle approximation would conclude that the "thin pancake" remains 
just outside the horizon. Our numerical simulation shows that the matter accumulates in 
a thin shell just inside the horizon, not outside. To observe this analytically would require 
one to include the back-reaction. The numerical simulations of spherically symmetric collapse 
previously mentioned [7]-[ll], do not observe the thin shell at the horizon discussed in |12j . 
This is not surprising since these studies do not use a coordinate corresponding to the time 
measured by a static observer at spatial infinity. 

The viewpoint of the static asymptotic observer is important for understanding black hole 
thermodynamics. The temperature T = 1/(8ttM) of a Schwarzschild black hole is the tem- 
perature as seen by an observer at rest at spatial infinity [12] and its entropy represents a 
measure of an external observer's ignorance of the internal configurations hidden behind the 
event horizon [13\ [6] . The formation of a thin shell near the horizon in isotropic coordinates 
plays a significant role in black hole thermodynamics and the holographic principle. In ther- 
modynamic studies [HE], it was observed that the black hole entropy is concentrated in a thin 
shell just inside the horizon. Isotropic coordinates offers a vantage point from where one can 
observe the entropy accumulate directly on the surface of the horizon and this represents a 
direct and dynamical observation of the holographic principle [14} [T5] . 

Our paper is organized as follows. In section 2 we state the equations of motion, initial states 
and boundary conditions that are used in our numerical simulation (these are derived in detail 
in [2]). In section 3 we discuss the formation of the black hole. In particular, we state the 
defining property of the N = spacelike two-surface and discuss how matter enters into the 
interior as the N = surface expands (and is identified as the apparent horizon at late times). 
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In section 4 we analyze the entire spacetime: the interior collapsing universe, the thin shell 
region and the exterior vacuum region, which acts as a benchmark for our numerical simulation. 

2 4D gravitational collapse of a scalar field 
in isotropic coordinates 

The spacetime during the collapse process is described by the spherically-symmetric time- 
dependent metric ([2]). The function ip(r,t), the conformal factor, is a dynamical variable 
but the lapse function N(r, t) is not 0- Though it is not dynamical, each time step yields a 
new spacelike hypersurface Sj so that one must reevaluate the value of N(r, t) at the new 
surface. This is governed by an ordinary differential equation instead of an evolution equation. 
The Einstein field equation governing the evolution of tp contains two time-derivatives. For 
numerical purposes, this is broken into two single time-derivative equations. The Einstein field 
equations also yield two constraint equations: one for the energy and the other for momentum. 
These are not evolution equations but are useful for monitoring the accuracy of the simulation. 
The collapsing matter is a Klein-Gordon massless scalar field x( r > t) which is minimally coupled 
to the gravitational field. Again, the evolution of x 1S governed by two single time-derivative 
equations. Once the equations of motion are supplemented with initial states and boundary 
conditions, the evolution of the metric and matter fields is unique. Depending on the initial 
state, the system either collapses to a black hole or disperses due to the internal pressure of 
the scalar field. We are interested only in the case of black hole formation. 

2.1 4D equations of motion 

The Einstein field equations yield an energy constraint, a momentum constraint, an evolution 
equation for the conformal factor ip, an evolution equation for its "momentum conjugate", the 
trace of the extrinsic curvature K and finally an ODE for the lapse function N . In 4D, the 

1 The lapse changes with time. By "not dynamical" we mean that its momentum conjugate is not defined 
because the gravitational action does not contain time-derivatives of N [4]. 
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equations are (see [2] for details and derivation): 

4 2 2 K 2 

energy constraint: — tp = k £ — . (4) 

K ' K 2 X , 

momentum constraint: — = y (5) 

3 2 N A v ; 

evolution equation for tp : — = — (6) 

K K 2 M t A , 3 2r/i N o f : 1 



evolution equation for K: — = -r-(l + w) + —K £(1 — w) — 6—j. ( — H — 

.«'^ + I)_ (ta ^* (7) 



iV?/; 4 V if) r ' Tp 

— f> ( ^' \ _L r f) ( N> \ - _ 2^' 2 



ODE for N: — d r ( + — d r ( — = -k 2 ^- (8) 

- ,2 r V ™/.3 / /V V i4 7 ,/,4 v 7 



where = 8tt G and the arbitrary parameter w in arises from a term added for numerical 
stability purposes (we set it to unity) . The energy density of the matter field is defined as 

There are also the evolution equations for the matter field x an d its "momentum conjugate" 
p. These are given by [2] 

p = ^ 6 § (10) 

and 



2.2 Initial states and boundary conditions 



We choose the same initial states and boundary conditions as in [2]. We will therefore be brief. 
We choose a static initial state, where K and p are zero. From © and (|10p . this implies that 
^ = x = at time t = 0. For this initial state, the momentum constraint ([5]) is automatically 
satisfied. The initial field configuration for the scalar field x is chosen to be [2] 

8A 2 r 4 

X(* = 0,r)= (A2 + r2)4 ) (12) 

where A is a scale parameter. Note that the though the scalar field configuration has one 
peak, the energy density profile has two peaks. For A sufficiently small, the self-gravitational 
attraction is large enough to initiate a gravitational collapse to a black hole. 
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The initial state for the conformal factor ip is obtained by solving the energy constraint @ 
with the static initial conditions (K =p = 0) i.e. 



2 



_ 5r (rV) = --xV (13) 



Asymptotically, the spacetime is flat: tp — > 1 and N — > 1. We therefore express ip as a power 
expansion in the parameter C = ft2 i- e - 

i> = 1 + C^i + CV2 + • • • • (14) 

and matching the left hand side with the right hand side of (|13p order by order in Q yields the 
initial state for ip [2]: 

ip = 1+A tt^ ) r^ g [l575r 15 A+12075r 13 A 3 +.. + 1575(r 2 +A 2 ) 8 tan- 1 (r/A) 



A 3 40320r (r 2 + A 2 ) 8 

+ 



r 20101331257r 2 - . „ 2 , ON1fl ( - r A 2 (2967339375 r 30 
r(r 2 + A 2 ) 16 



A« 10538886758400 
+47828405625 r 28 A 2 + . . . + 60640845881 r 2 A 28 + 41 13561991 A 30 ) 
+1276275 A(r 2 +A 2 ) 8 (825 r 16 +4500 r 14 A 2 + . . + 15436 r 2 A 14 + 2717A 16 ) tan _1 (r/A) 

+4020266250r (r 2 + A 2 ) 16 [tan~ 1 (r/A)] 2 ) 
+ ... 

Our expansion is carried to six orders but we only show here the first two order terms. Once 
the initial state of ip is obtained, the initial state for iV is obtained numerically by solving its 
associated ODE (jHJ). We iterate backwards to obtain iV starting at the outer computational 
boundary r = R, where the spacetime is flat and N = l. 

The above equations and initial states need to be supplemented with appropriate boundary 
conditions at the origin r = and the outer computational boundary r = R. To ensure 
regularity of the solution at r = 0, we impose the following boundary conditions: x'(0,t) =0 
and i^'(0,t) = 0. To ensure asymptotic flatness at the computational outer boundary r = Rwe 
impose N(R, t) = 1, K'(R, t) = and p'(R, t) = 0. Together with the initial matter and metric 
states, these lead to a unique evolution. 



3 Black hole formation 

We work in geometrized units where G = c=l (energy, mass, time and distance are measured 
in length). 
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As in [3j, a black hole will be considered to have formed when the lapse function N(r, t) crosses 
zero. Let ro be the radius where N = 0. If this radius remains fixed with time, light could 
not cross the N = two-surface in either direction since it takes an infinite amount of time t 
for light to travel from any point to this fixed surface (as viewed of course by an asymptotic 
observer whose clock measures the time t). However, in our collapsing system, the radius ro 
changes as a function of time and increases until it reaches a radius of r^ = 0.335 at late times 
(this is the radius in figure [T]) where the lapse function crosses zero at t = 22). Note that N 
remains negative inside once it crosses zero. The defining property of the N = two-surface 
during the simulation is that outgoing radial null geodesies inside the surface never cross it 
as it evolves. The radius as a function of time for outgoing null geodesies that start inside 
at different radii is plotted in figure [2j The radius tq as a function of time is plotted as well. 
The null curves never intersect the ro curve; they never cross the N = surface as it evolves. 
The null curves become flat at late times and converge towards ro reaching a radius just below 
rh = 0.335. Outgoing (and of course ingoing) null geodesies inside are trapped within this 
radius. 

A simple analytical argument also shows that outgoing null geodesies on the inside do not cross 
the N = surface. From the metric, radial null geodesies obey the relation [dr/dt) 2 = N 2 /^ 
and this implies that dr/dt = at the N = surface (for both ingoing and outgoing null 
geodesies). The condition for radial null geodesies to cross the surface from inside to outside 
is that dr/dt — dro(t)/dt be positive when evaluated at r = ro- This condition is clearly not 
satisfied since dr /dt = at r = ro and dro(t)/dt > 0. However, an ingoing radial null geodesic 
can cross the surface from outside to inside. This requires that dr/dt — dr^/dt be negative 
when evaluated at r = ro and this is automatically satisfied since again dr/dt = at r = ro 
and dro/dt > 0. In particular, matter outside the surface can penetrate inside allowing the 
black hole mass to increase during the evolution. At late times, the N = surface is basically 
stationary at r = ry l (drq/dt — > 0) and the flow of matter from outside to inside comes to a halt 
with the black hole mass reaching its maximum value of Mbh- 

The expansion scalar G for outgoing null geodesies orthogonal to the spacelike two-surface 
N(r,t) = (the two-sphere ds 2 = ^ 4 r 2 df2 2 ) can be readily calculated [16] and is equal to 



where = g^ u + + £^k u is the two-dimensional metric and fc M = — — (l/|iV| , l/?p 2 ,0, 0) 
1 

and = —={1/\N\, — l/^ 2 , 0, 0) are future-directed null vectors orthogonal to the two-surface 
V2 

with normalization k^i^ = — 1. is tracked numerically, and at late times, it is equal to zero at 
the iV = surface. In particular, the spacelike two-surface r = r^ = 0.335 is an apparent horizon. 




(15) 
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The event horizon is a null hypersurface which is defined globally and requires knowledge of 
the entire future history of the spacetime. In practice, however, it can be located after a 
finite evolution time. During the simulation, we plotted the norm of the vector orthogonal to 
hypersurfaces of constant N. We identified a null surface at a radius r nu u outside the N = 
surface {r nu u > r ). With time, r nu u and ro approach each other so that the event horizon 
is located at a radius just outside the apparent horizon (r/j = 0.335). This is in agreement 
with the fact that outgoing (and of course ingoing) null geodesies inside the event horizon are 
trapped and cannot escape out to infinity (see fig. ([21)). 

We know from the exterior analytical solution in isotropic coordinates that = Mbh/ 2 where 
Mbh is the black hole mass. In our dynamical collapse scenario, the black hole mass Mbh 
is less than the total (conserved) ADM mass because part of the ADM mass stems from an 
outgoing matter wave. The function m(r, t) = —r 2 tp' evaluated at infinity corresponds to the 
ADM mass [HE]- To illustrate the difference between the black hole mass and the ADM mass, 
we plot m(r,t) as a function of r at different times in figure [3j Note that at large r, one can 
distinguish between two plateaus at late times. The value at the first plateau corresponds to 
the black hole mass Mbh and the value at the second plateau is the ADM mass. The difference 
between the two values is due to an outgoing matter wave which propagates to infinity with 
time. The identification of the first plateau with Mbh is confirmed by the fact that Mbh/2 
at late times is equal to = 0.335 to within 2%. In figure El the function m(r, t) at late 
times dips to a negative value inside the horizon and then rises to reach a final positive value 
outside the horizon. The origin of these two contributions can be determined by expressing 
the function m(r,t) in integral form via the energy constraint (jj]): 



The negative dip stems from the — K 2 term, a gravitational term, and the rise to a positive 
value stems from the positive energy density £ of the matter. Their sum (the integral) yields 
the positive mass Mbh of the black hole. This is a clear illustration of the fact that the mass 
of a black hole stems not only from matter but also from gravitation itself jlj. 

4 Features of the spacetime: numerical results 

4.1 Exterior vacuum region: a benchmark 

Past numerical studies [H [2] have shown that the collapse of matter in isotropic coordinates 
leads to metric and matter functions whose gradients change sharply in the region near the 
apparent horizon. To place more points in the interior region we expressed the radius r as a 
function of a parameter x: r = x runs from to 1 and r runs from to "infinity" (i.e. the 




(16) 
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outer computational boundary). At highest resolution, we used 2 x 10 points for a step size of 
Ax = 5 x 10~ 5 (and At = 1 x 10" 4 ). We used the initial state of {T2]) with A = 1.5. To evolve 
the fields in time, we used a fourth-order Adams-Bashforth-Moulton (ABM) explicit scheme. 
During the evolution we monitored the ADM mass i.e. M = 4tt/ c f{£ - l-r)^ 5 r 2 dr where R 
is the outer computational boundary, £ is the energy density given by © and K is the trace 
of the extrinsic curvature (k 2 = 8ir since we work in geometrized units where G = c=l). 

The metric for the static Schwarzschild vacuum exterior in isotropic coordinates is known 
analytically and is given by ([3]). The numerical simulation should reproduce this metric in the 
exterior region and this provides a test on the accuracy of our code. In particular, the conformal 
factor in the exterior should have a form close to ip = 1 + Mj (2r) and its peak value should 
be close to 2. The theoretical and numerical curves for ip are shown in figure HI Note how 
the numerical curves approach the analytical curve with time. At late times, the percentage 
difference between the two curves at any point is always under 1.31%, being significantly under 
this value most of the time. The peak value of ip analytically is equal to the integer 2. For 
the run mentioned above we obtain a peak numerical value of 1.974, when the ADM mass 
deviated by 1% (at t = 18.28) and a peak value of tp of 1.992 when the ADM mass deviates 
by 5% (at a later time t = 20.26). This corresponds to a tiny difference of 1.3% and 0.4% 
respectively with the expected analytical value. Note that the code is still evolving well (still 
approaching the analytical curve) from t = 18.28 to t = 20.26 even though the ADM mass 
has started to deviate from its original value. We ran the code for the following grid sizes: 
Ax = {10~ 3 , 5 x 10~ 4 , 2 x 10~ 4 , 1 x 10~ 4 , 5 x 10~ 5 }. Let fa and fa be the peak values of V at 
ADM mass deviation of 1% and 5% respectively. Their values for the corresponding grid sizes 
are fa = {1.783, 1.859, 1.923, 1.953, 1.974} and fa = {1.897, 1.934, 1.966, 1.982, 1.992}. In both 
cases the values are converging and approaching closer to the expected analytical value of 2. 

We also ran the code with different initial states: different values of A for the initial state of 
the form (|12p and a different initial state of "non-Gaussian" form. The results are not very 
different and the features of all graphs remain the same. At highest resolution, we obtained 
fa = 1.970 and fa = 1-989 for A = 1.55 and fa = 1.966 and fa = 1.986 for A = 1.60. For a 
"non-Gaussian" initial state of the form \ = 1+ "^ 1 2r 4 with Ai = 0.3 amd A2 = 0.05 we obtained 
fa = 1.972 and fa = 1.978, again at highest resolution. 



4.2 Thin shell structure 

At late times, the energy density is concentrated near two peak^l just inside the apparent 
horizon, the iV = spacelike two-surface at r/j = 0.335 (see fig. [5]). During the collapse, the two 
peaks approach each other and the energy density profile becomes thinner with time so that 

2 Recall that the initial scalar field configuration (|12fl has an initial energy density profile with two peaks. 
The two peaks remain throughout the evolution. 
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the two effectively merge into one thin shell located just inside the horizon. The trace of the 
extrinsic curvature, K, takes the form of a Heaviside function as shown in fig. [61 being very 
close to zero outside the horizon and being negative and spatially constant inside. The jump 
discontinuity in the extrinsic curvature is precisely what is expected from the presence of a thin 
shell of matter at the location where the jump occurs [3J. An ideal thin shell is expected to have 
a delta- function singularity in the stress-energy tensor and the Riemann tensor Ry ar [I]- 
In our numerical simulation, these singularities are represented by various quantities increasing 
with time and reaching peak values or spikes at late times. We have already looked at the 
energy density (fig. [5]), which is basically zero everywhere except for spikes located just inside 
the horizon and whose peak value increases with time. The Ricci scalar behaves in roughly the 
same fashion: it is zero everywhere except for spikes which increase in magnitude with time in 
the thin shell region (see fig. [7|). The square of the Weyl tensor (C 2 = C^ aT C^ uar ) has a few 
peaks just inside the horizon which can be seen to increase significantly with time (see fig. 
The spacetime can therefore be said to be evolving towards a curvature singularity at the thin 
shell (note that this occurs as the shell "falls" towards the origin. See paragraph below). There 
is both a Ricci singularity associated with a spike in the energy density (and the Ricci scalar) 
and a Weyl singularity associated with a spike in C 2 representing strong tidal forces. Note that 
the standard Schwarzschild spacetime has a Weyl singularity but no Ricci singularity, since it 
is a purely vacuum solution. In the prototypical Oppenheimer-Snyder collapse [T7], modeled 
after a contracting FLRW universe, the curvature singularity is due to a high energy density; 
there is a Ricci singularity but no Weyl singularity since the Weyl tensor is identically zero for 
an FLRW universe. Our collapse scenario leads to a more physically realistic situation where 
there is both a Ricci and Weyl curvature singularity. 

We will see shortly that the interior is a collapsing universe where ip is almost homogeneous 
and approaches zero. The thin shell is therefore "falling" towards or approaching r = while 
remaining at a radial coordinate near r^. The proper radial distance between r = and r = 
is given by Jq H ip 2 (r,t)dr which is approximately equal to "0 2 (t)r/j (since ip is approximately 
homogeneous in the interior). This tends towards zero as tp approaches in the interior. The 
curvature scalars evolve towards the singularity as the shell "falls" towards the origin. 

4.3 Interior spacetime: a collapsing isotropic vacuum universe 

Besides the emergence of a thin shell structure, the other salient feature of the spacetime is the 
interior. The interior is defined here as the nearly zero energy density (almost vacuum) region 
inside the apparent horizon where the thin shell of matter is absent. In the interior ip and N 
are close to being homogeneous (see figures Q] and |4| so that ip(r,t) ~ ip(t) and N(r,t) ~ N(t). 
The magnitude of ip and N are also decreasing with time towards zero: N is approaching 
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zero from below and if) from abov^l. The interior can be viewed approximately as a collapsing 
FLRW universe (approximately because the interior is not perfectly homogeneous but has a 
slight inhomogeneity) . Defining r = — j N(t)dt the metric can be cast approximately in the 
FLRW form with flat spatial sections: ds 2 ~ —dT 2 + i[) 4: (T)(dr 2 + r 2 d0 2 jj|. The point of casting 
the metric in this approximate form is to extract the Hubble parameter as a measure of the 
rate of collapse. The scale factor is a(r) = iP 2 (t) and the Hubble parameter is H = = 

= — -^jf where a dot represents derivative with respect to the original t coordinate. The 
Hubble parameter is plotted as a function of time in figure [9j It increases rapidly in magnitude 
with time and is negative since the interior is collapsing or contracting. The interior is clearly 
heading towards a curvature singularity as ip — > 0. The areal radius tp 2 r approaches zero inside 
and the Weyl tensor squared (C 2 ) in the interior region increases rapidly with time (see fig. 
EH. 



The trace of the extrinsic curvature K = —6ip/(Nijj) can be expressed in terms of the above 
Hubble parameter as K = 3H. The energy constraint equation @ then takes the form 

In a perfectly homogeneous spacetime as in FLRW, V 2 iJj would be zero and the above would 
reduce to the Friedmann equation for flat spatial sections and zero cosmological constant: 
H 2 = ^f-£ where £ is the energy density of the matter. Our interior is not exactly FLRW: 
it has a slight inhomogeneity which implies that \7 2 t/j ^ 0. Also, the interior tends towards a 
vacuum: £ — > 0. For the interior, the governing equation is therefore H 2 = ^5 \7 2 ip. We have 
seen that the magnitude of the Hubble parameter increases rapidly with time. Though there 
is only a slight inhomogeneity in ip so that V 2 ip is small, ij) tends to zero in the interior so that 
the right hand side (RHS) term, which contains tp 5 in the denominator, increases rapidly with 
time also. The RHS term also has a negligible spatial dependence like K and H. 

In the exact FLRW case governed by the above Friedmann equation, the energy density even- 
tually diverges as in Oppenheimer-Snyder collapse [T7]. In our slightly inhomogeneous con- 
tracting interior vacuum spacetime, it is the term V 2 ^ that eventually diverges as tp — > 0. 
The slight inhomogeneity plays a crucial role in the collapse process. The curvature singularity 
in the exact FLRW case is a Ricci singularity (due to high energy density). In contrast, our 
collapsing interior vacuum spacetime can be said to be evolving towards a Weyl curvature 
singularity, the same type of singularity encountered in the Schwarzschild vacuum interior. 

3 In the interior, iV and i/j evolve towards zero but are never identically equal to zero. 
4 The collapse takes an infinite time t but only a finite proper time t. 
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5 Conclusion 



We studied numerically the spacetime which emerges after the collapse of a massless scalar 
field in isotropic coordinates. A salient feature of the spacetime is the emergence of a thin shell 
of matter just inside the horizon. The spacetime consists of two vacuum regions separated 
by the thin shell with the Ricci scalar being zero everywhere except for a spike at the shell. 
There is a jump discontinuity in the extrinsic curvature across the horizon, in accord with the 
presence of a thin shell near the horizon. The expected delta function singularities at an ideal 
thin shell [1] are represented numerically as spikes in the Ricci scalar, the energy density and 
the Weyl tensor squared which all increase significantly with time in the thin shell region. The 
exterior region settles into the static Schwarzschild vacuum solution in isotropic coordinates. 
The interior is an isotropic universe that collapses basically in vacuo. The magnitude of the 
"Hubble parameter" and the trace of the extrinsic curvature increase with time not because 
the energy density increases, but because a term proportional to V 2 ^/^ 5 , which has a non-zero 
value due to a slight inhomogeneity, grows as tp tends to zero. The Weyl tensor squared rises 
rapidly with time in the interior as the areal radius shrinks. 

The most natural case to study next would be charged collapse in isotropic coordinates. There 
are many interesting questions to explore. Would the charge density be concentrated in a thin 
shell and have a similar profile as the energy density? Would there be a timelike singularity 
and inner horizon as in the Reissner-Nordstrom (RN) spacetime? This question is important 
because the inner horizon in RN is thought to be unstable [1] and charged collapse studies 
in other coordinates have not revealed a timelike singularity |10^ 111]. But more importantly, 
isotropic coordinates would allow one to study the thermodynamics of the charged black hole 
at late stages of the collapse. In particular, one could track the gravitational Lagrangian, a 
function which has been shown to be connected to thermodynamic properties of a black hole 
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r 



Figure 1: Lapse function N as a function of r at different times. The radius ro where N = 
increases with time and reaches = 0.335 at late times. The two-surface at is an apparent 
horizon. Note that iV is negative in the region r < tq and approaches zero from below. 
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Figure 2: The radius of outgoing null geodesies as a function of time. The solid line represents 
ro. ri stands for the position of the photon particle at t = 10. The particles start inside the 
iV = surface where ro = 0.30 at t = 10. The null curves never cross the ro curve during the 
evolution; they remain inside the N = surface as both evolve. 
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Figure 3: The mass accumulation function m(r, t) as a function of r at different times. The 
first plateau at m ~ 0.67 corresponds to the black hole mass Mbh and the second plateau 
corresponds to the ADM mass. 
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Figure 4: The conformal factor ift as a function of r at different times. The solid line represents 
the exterior Schwarzschild analytical solution. 
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Figure 5: The energy density £ as a function of r at different times. 
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Figure 6: The extrinsic curvature K as a function of r at different times. The jump disconti- 
nuity is due to the presence of a thin shell of matter near = 0.335. 
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Figure 7: The Ricci scalar R as a function of r at different times. At late times R is basically 
zero everywhere except for a spike near the apparent horizon. 
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Figure 9: The Hubble parameter as a function of time. The data is taken at r = but H has 
almost no dependence on r. 
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Figure 10: The Weyl tensor squared C 2 as a function of r at different times with a focus on the 
interior region (smaller r values). Its value increases significantly in a very short time interval. 
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